Relativistic quantum theory of high harmonic generation on atoms/ions by strong 

laser fields 
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High-order harmonic generation (HHG) by hydrogenlike atoms/ions in the uniform periodic elec- 
tric field, formed by the two hnearly polarized counterpropagating laser beams of relativistic inten- 
sities, is studied. The relativistic quantum theory of HHG in such fields, at which the impeding 
factor of relativistic magnetic drift of a strong wave can be eliminated, is presented arising from 
the Dirac equation. Specifically, a scheme of HHG in underdense plasma with the copropagating 
ultraintense laser and fast ion beams is proposed. 

PACS numbers: 42.65.Ky, 32.80.Rm, 31.30.J-, 03.65.Pm 



I. INTRODUCTION 

The unprecedented progress of laser technology in the 
recent decade made available the realization of fem- 
tosecond laser sources vifith relativistic: 10^* W/cm'^ 
(in optical domain) and ultrarelativistic intensities up to 
10^2 \Y/cm^ by the chirped pulse amplification technique 
It is expected that the next generation of ultrarel- 
ativistic lasers will significantly exceed these intensities 
in the near future [2| . The interaction of such powerful 
electromagnetic radiation with matter opens real possi- 
bilities for revelation of ultrafast laser-matter interaction 
dynamics in supershort time scales with new research 
fields, such as Relativistic Optics 0-01 ^"^^ Attoscience 
@-8i], as well as many nonlinear electrodynamic efi^ects 
at the ultrarelativistic interaction with the QED vacuum 
and accelerator beams (9||. 

At the laser-matter interaction the high-order har- 
monic generation (HHG) by atoms/ions pi]| - [l3j is a sub- 
ject of extreme interest, since HHG is one of the most 
promisiiig mechanisms for creation of coherent XUV 
sources [H, [l^ and intense attosecond pulses Q , in par- 
ticular, for the control of electrodynamic processes on 
attosecond time scales. To reach a far X-ray region one 
needs the atoms or ions with a large nuclear charge and 
laser fields of ultrahigh intensities at which the nondi pol e 
interaction and relativistic effects become essential |16l - 
[Tgj . At such intensities the state of ionized electron be- 
comes relativistic already at the distances Z ^ A (A is 
the wavelength of a laser radiation) irrespective of its 
initial state, hence the investigation of laser-atom/ion 
induced processes such as Stimulated Bremsstrahlung 
(SB), Above Threshold Ionization (ATI) etc., require rel- 
ativistic consideration [20l - l22j . Specifically, the relativis- 
tic drift of a photoelectron due to the magnetic field of a 
strong electromagnetic wave becomes the major inhibit- 
ing factor in the relativistic regime of HHG. Due to this 
drift the significant HHG suppression takes place and di- 
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verse schemes have been proposed to overcome this neg- 
ative effect [l^-dHl- For compensation of the magnetic 
drift in relativistic domain, in Refs. [1^, [l^l two coun- 
terpropagating laser pulses of linear and circular polariza- 
tions were considered, respectively, at which the configu- 
ration of a standing wave is formed. The effect of result- 
ing magnetic field of the standing wave is vanished near 
the stationary maxima at the waves' linear polarizations 
(23} . While at the circular polarizations of counterprop- 
agating laser pulses one can achieve the fully vanishing 
of resulting longitudinal magnetic force, responsible for 
magnetic drift [24!|. However, the latter takes place at 
the adiabatic turn on of the waves [9|] , which is not valid 
for short laser pulses of relativistic intensities. 

So, a standing wave configuration formed by the two 
counterpropagating laser beams of linear polarizations is 
of interest due to the simplicity to realize such field struc- 
ture providing incomparable large HHG rates in the rela- 
tivistic regime. At the lengths much smaller than a wave- 
length of a pump wave the effective field of the standing 
wave may be approximated by the uniform periodic elec- 
tric field. Note that the consideration of the HHG in the 
field of a standing wave in the paper [l^ is based on the 
semiclassical model, and only the first order relativistic 
effects have been taken into account. The entire relativis- 
tic quantum-mechanical consideration of the HHG in a 
standing wave or in a nonstationary strong electric field 
on the base of the Dirac equation is still lacking. 

To overcome the suppression of the HHG in the rela- 
tivistic regime caused by the magnetic field of a strong 
laser pulse, in the present work the relativistic HHG by 
hydrogenlike ions in the uniform periodic electric field, 
formed by the two linearly polarized counterpropagat- 
ing laser beams of relativistic intensities, at the distances 
much less than a laser wavelength, is studied. We pro- 
pose also a scheme of relativistic HHG in plasma, where a 
traveling laser pulse in the own frame of reference (mov- 
ing with the pulse group velocity in plasma) is trans- 
formed into the uniform periodic electric field, i.e. the 
wave magnetic field in this case vanishes completely, in 
contrast to a standing wave configuration formed by the 
counterpropagating laser beams. So, it is also feasible 
the effective HHG in plasma with the copropagating ul- 
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traintense laser and fast ion beams of the same velocities. 
Note that in Refs. [1^ H^] a scheme with counterprop- 
agating relativistic ion and laser beams was considered, 
where due to the large Doppler upshift the high harmonic 
frequency may appear in the 10 keV-1 MeV domain. 

The organization of the paper is as follows. In Sec. II 
we present analytic and fully relativistic quantum theory 
of HHG arising from the Dirac equation. In Sec. Ill with 
the help of Fast Fourier Transform algorithm we present 
HHG spectrum for diverse ion and field parameters. Fi- 
nally, conclusions are given in Sec. IV. 



II. BASIC MODEL AND THEORY 

Let two linearly polarized plane electromagnetic waves 
with carrier frequency w and amplitude of the electric 
field Ea 

El = Ea cos(a;t — kr), E2 = Ej, cos(cj< + kr), (1) 

propagating in the opposite directions in vacuum, inter- 
play with the hydrogenlike ions having the charge number 
of the nucleus Za- We will assume that A >> a, where a 
is the characteristic size of the atomic system and A is the 
wavelength of a pump wave (for the HHG this condition 
is always satisfied). 

At the photoionization of an atom/ion in the strong 
traveling wave field taking into account the relativistic 
drift due to the magnetic field, one can expect that the 
probability of ionized electron recombination with ionic 
core could be non-negligible only if the electron is ini- 
tially born with a nonzero velocity oppositely directed to 
incident laser beam. The probability of tunneling ion- 
ization with nonzero initial velocity is given by quantum 
mechanical tunneling theory [2^ as: 

t^,o„oce-^(^"+^^r/(^^), (2) 

where v is the initial velocity, and E is the electric field 
strength of the wave (here and below, unless stated oth- 
erwise, we employ atomic units). However, according to 
([2]) the ionization probability falls off exponentially if this 
velocity v becomes larger than the characteristic atomic 
velocities, irrespective of it's direction. Since we study 
the case of superstrong laser fields with ^ = S/ca; ~ 1 
(^ is the relativistic invariant parameter of the wave in- 
tensity) when the energy of the interaction of an electron 
with the field over a wavelength becomes comparable to 
the electron rest energy, the required velocities becomes 
comparable with the light speed c (c = 137 a.u.). Hence, 
the probability ^ in the such fields is practically zero. 
Therefore, in considering case of a standing wave formed 
by the laser beams H]), a significant input in the HHG 
process will be conditioned by the ions situated near the 
stationary maxima of the standing wave. For this points 
the magnetic fields of the counterpropagating waves can- 
cel each other. Since the HHG is essentially produced at 



the lengths Z << A, we will assume the effective field to 
be: 

E(t) = e£;ocosa;i, (3) 

where = 2Ea, and e is the unit polarization vector. 

As is clear from the above consideration for HHG, 
one needs to exclude the relativistic drift of a photo- 
electron due to the magnetic field of a traveling wave, 
i.e. the magnetic component of the wave. This can also 
be achieved in the plasma- like medium with a refractive 
index np{uj) < 1. Indeed, let a plane, transverse, and 
linearly polarized electromagnetic wave with a frequency 
uj propagates in plasma {uj > ujp; u!p — ^AnNg is the 
plasma frequency). For the electric and magnetic field 
strengths we have 

E = Ea cos {ut -kr), H =- [k x E] , (4) 

where |k| — np{ijj)Lo / c. Suppose that an ion beam co- 
propagates together a laser beam in such plasma with 
a mean velocity V equal to a laser beam group velocity: 
V = cnp{(jj). To find out the HHG rate it is convenient to 
solve the problem in the center-of-mass (C) frame of ions. 
In this frame the wave vector of the photons is k' = 0. 
The traveling electromagnetic wave is transformed in the 
C frame into the field [9]: 

H' = 0, E'(t') = E^jCosw't', 



Uj' = u^l - n'^{uj) = UJp. (5) 

Hence, in the frame of ions center-of-mass the prob- 
lem of HHG in a plasma is reduced to one in vac- 
uum in the field of the same configuration (jSj with 
Ejj = EaWp/w. When oj » ujp, one needs relativistic 
ion beams. Concerning the ion beams of necessary rela- 

tivism 7 = (1 —V^/c^) — uj/ujp, note that the rel- 
ativistic ion beams in arbitrary charge states, with the 
Lorentz factor up to about 30 is supposed to be real- 
ized at the new accelerator complex at Gesellschaft fiir 
Schwerionenforschung (GSI) (Darmstadt, Germany). It is 
worthy to note that if a laser frequency w is close enough 
to plasma frequency Wp, one can achieve the implemen- 
tation of proposed scheme for HHG with nonrelativistic 
ion beams. In this case one can use quasimonoenergetic 
and low emittance ions bunches of solid densities gener- 
ated from ultrathin foils -nanotargets by supershort laser 
pulses of relativistic intensities [29|. 

To find out the relativistic probabilities of HHG in the 
field ([3]) we arise from the Dirac equation: 

^^^(^cap + ^c' + ^-rE{t)y^), (6) 

where a and (3 are the Dirac matrices in the standard rep- 
resentation, a = i(Jx,cFy,az) are the Pauli matrices and 
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p is the operator of the kinetic momentum (p = — iV). 
Without loss of generahty one can take the polarization 
vector e aligned with the z axis of spherical coordinates. 

We denote the atomic bound states by [yy), where 77 
indicates the set of quantum numbers that character- 
izes the state 77 = Here n is the principal 
quantum number, j is the whole moment, I is the orbital 
moment and M is the magnetic quantum number. 

Following the ansatz developed in the Ref. [l^, the 
time dependent wave function can be expanded as 



(7) 



Here F (x) is the Euler gamma function, 6 and ip are 
the polar and azimutal angles, e = 1 — a/1 — Z^/c^, and 
e — (1 — e) is the energy of the ground state. We 
assume Za < 137, and the parameter e can take values 
lying in the interval < e < 1. In the expansion ([7]) 
we have excluded the negative energy states, since the 
input of the particle-antiparticle intermediate states will 
lead only to small corrections to the processes considered. 
Neglecting the depletion of the ground state Co {t) ~ 1 
and the free-free transitions due to Coulomb field, the 
Dirac equation for {p,t) reads as: 



dC^ {p,t) 
dt 



E{t) (p)-e)C^(p,t) 



Here 



|P,Ai) 



1 f(p)+c2 



(27r 



,3/2 ' 



2£ (p) 



(8) 



are the Dirac free solutions [3^ with energy £ (p) — 
\/c?p^ + c? and polarization states /i = 1, — 1: 

^1 = {]^ . ^-1 = f?) ■ (9) 



As an initial bound state wave function I^q) we assume 
the ground-state bispinor wave function for the hy- 
drogenlike ion with the quantum numbers n — 1, j — \ /2, 
I = 0, and M = 1/2: 



1%; = 



zl'' 



v/^r (3 - 2e) 



/ \/2^ 


zcos0^/e 
sin 6'e"^yiy 



(10) 



iV^ {p)E{t) + 



c^E jt) (/ip^ - ipy) 
2£ (p) {£ (p)+c2) 



C_^(p,t), (11) 



where 



V^{p) = {p,fi\z\7jo) (12) 
is the atomic dipole matrix element for the bound-free 
transition. The latter can be calculated with the help of 
integral [Slj : 



ici cos 9 cos 



Jm (ci sin e sin e)P/" (cos 9) sin eide* 



1/2 



— ) 7'-"Pr(cOSe)J,+ i/2(ci), 
Ci ' 



where P™ (cos O) is an associated Legendre function of 
degree / and order m, Jm is a Bessel function of order m. 
Thus, for (p) we have: 



Pi (p) = i- 



23-^ /^T^Z^/'-'F (3 - e) 



2f v/r(3-26) (p2 + Z2)3- 



1 ^g'-^ 

4cv^ 2p2 



^3 — e '^13- 



1 iZ, 



3-e 



(2 - e) 



f + c2 Y 2 (2 - e) 



_ "ft3-eN _ 2» (p^ + Z^) / 2-c I "(-t2-e^ _ 2F (1 e) (p^ + ^g) / _ 

^ ^ (2-e)Zap ^ ^ ^ r(3-e) Z.V ^ 



(13) 



(P) = 



21+^ V v/r (3-2e)£:(£: + c2) 



^ ^"'^ j Cy^-'^ ytl-e') ^apF (2 ~ e) | -^12- 



(p2 + Z2)^ 



(p2 + Z2)^ 



, (14) 



where 



The last term in Eq. (|lip describes spin fiip due to 



T = 1 -i^, = 1-l-i 



Za 



(15) 
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free-free transitions and proportional to transverse mo- 
mentum Px,y As will be seen below, the main contribu- 
tion to the HHG process is conditioned by the electrons 
with momentum components p^^y — 0. Hence, we will 
neglect the last term in Eq. (fTTj) . which is valid for the 
fields E « Es, where Es — rin?c' /eh is the Schwinger 
field strength. Then, from Eq. pT|) for the probability 
amplitudes (p,0 we obtain: 



In contrast to nonrelativistic [T^ and relativistic [T3| 
cases for a traveling wave, this equation can not be solved 
analytically and requires numerical solution. Integrating 
the latter over p, for the current density we obtain: 



J{t) = {2nf'^Vl / dt'- 
Jo 



-i(S(pa,t,t')-e(t-i')) 



C, (P,0 = ^ dt'V, (p + i (A (t) -A (O)) E (f) ^ (^p^ _ 1^ 



X exp 



£ p+-{A{t)-A{t"))]-e 



dt" 



(16) 

where A {t) = — ec^^ sin ijjt is the vector potential of re- 
sulting electric field. 

For the harmonic radiation perpendicular to the polar- 
ization direction e one needs mean value of the z compo- 
nent of the electron current density J {t) — c 3^ j^*). 
Using Eqs. (O, and neglecting the contribution by 
free-free transitions, we obtain: 



dp I dt'J^ 





where 



xV^[p--A{t')]E{t') 



exp {-iS (p, t, t') +ie{t- t')} + c.c, (17) 



S{p,t, t') ^ j\ (yP-\^ (t")^ dt" 

= ^ {p+tc^smujt"f + c^dt" (18) 

is the relativistic classical action of an electron in the 
field and Jf, (p) =0(770! \p,fJ.)- 

As in the nonrelativistic case, the HHG rate is mainly 
determined by the exponential in the integrand of Eq. 
(IT7)) with exact relativistic classical action. The integral 
over the intermediate momentum p and time t' can be 
calculated using saddle-point method. The saddle mo- 
mentum is determined by the equation: 



dS{p,t,t') 
dp 



0, 



(19) 



which for momentum components gives p^;,!, — 0, and the 
z component of momentum (ps) is given by the solution 
of the equation: 



/' 

Jt' 



Ps+c^ sin Lot" 



-.dt" = 0. 



*' -J {ps+cS^ sin cut'' 



(20) 



where 



det^^'p 



dr 



-E{t') 



+ C.C., (21) 



dr 



Jt' lit',T-Oj Jt' 7Ht',T;0' 
7(i',r;e) = {l + {ps+c^sinujTf /c''y^\ (22) 

Note that contribution from spin flip due to bound-free 



transitions in the saddle-point vanishes: 2?- 



0. 



At ^ << 1, for detS'pp one can recover nonrelativistic 

resuh: det S'^'p = (< - t'f. 

The complex saddle times tg are the solutions of the 
following equation: 



dS{ps,t,t') 
dt' 



e = 0, 



(23) 



which may be expressed by the transcendental equation 



c^iPs+c^sinujtsf +c'^-e = 0. (24) 

Then expressing the saddle time as tg — th + i5, with 
lj5 << 1, one can obtain the saddle momentum 

Ps — —cS^sinLotj), (25) 

and the imaginary part of the saddle time: 



\E[h)y 



(26) 



where E (ti,) = c^cjcoswtf,. Taking into account Eq. 
([25|l . from Eq. (|20p for the real part of the saddle time 
we obtain 



sinujt" — sin ujti) 



-.dt" = 0. (27) 



1 + 1^2 {sinujt" — siniotbY 



As usual, if, is interpreted as the birth time of the pho- 
toelectron which returns at the moment t to the core 
and generates harmonic radiation. The transition dipole 
moment has singularity at the saddle times and the inte- 
gration has been made with the help of the formula: 

g{x)- -^dx ~ I'^gixo) [2/" (xq) A] = 



(X - Xq) 

r(iy/2)e-^^(^o) 



(28) 
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FIG. 1: The photoelectron energy gain in units of pondero- 
motive energy Up = and ionization time versus the 



Thus, taking into account Eqs. (|19II28I) . we obtain the 
uhimate formula for current density: 

= H ^ion ih) Cpr {t, h) Crec {t, U) + C.C. (29) 
tb 

The formula ((29)) is the analogous to the nonrelativistic 
formula for the dipole moment in the three step model 
[13I I . Here the summation is carried out over the solutions 
ofEq. (1271) . The tunneling ionization amplitude Cion (ife) 



C'ion {tb) 



2, r(^) \E{U)\i 
Vr(3 -2e)2*S(ifc) 



3-3e 



(30) 



The propagation amplitude is given by the expression 
(27r)^/^ exp{-iS' {p,,t, h) + ie{t- 4)} 



Cpr (i, tb) 



(31) 

electron time evolution in the continuum (return time) for and the recombination amplitude is: Crcc{t,tb) = 
the various laser intensities. J'^ (p^ _ 1\ (t)), where 



^'^P^ 2£ V2r(3-2e)(p2 + z2)2- 



(f +C2)p 




. r(l-e)(p^ + Zg) 
1 - e pZa 



(32) 
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III. RELATIVISTIC HARMONIC SPECTRA 
FOR DIVERSE ION AND FIELD PARAMETERS 

In this section we present some numerical simulations 
for HHG when relativistic effects become essential. In 
considering relativistic case the saddle time and energy 
gain essentially depend on the intensity of the pump 
wave. The latter leads to the modification of the HHG 
spectrum compared with nonrelativistic one. In Fig. 1 
we present the solution of Eq. (|27p for the born times tb 
which are limited to a quarter of the laser period. The 
Fig. 1 also illustrated the photoelectron energy gain ver- 
sus the electron's time evolution in the continuum (return 



time) for the various laser intensities. It is well known 
that the cutoff frequency of HHG is defined by the max- 
imal classical energy gain. For this reason in Fig. 2 the 
maximum energy gain of the photoelectron as a function 
of the relativistic parameter of the wave intensity is 
displayed. As we see, the relativistic cutoff essentially 
differs from the nonrelativistic one for ^ > 1 and the 
shift of the cutoff position to the lower values of the har- 
monic order for the same laser intensity becomes evident. 
The Fig. 1 also reveals the multiplateau character of the 
harmonic spectrum like to the nonrelativistic one. 

To find out harmonic spectrum, the Fast Fourier 
Transform algorithm has been used. 
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FIG. 2: Maximum energy gain of photoelectron, which defines 
the cutoff frequency as a function of the relativistic parameter 
of the wave intensity. 
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FIG. 3: Harmonic emission rate via logj^Q (|J]v|'^), as a func- 
tion of the harmonic order, for an ion with Za = -i, £, = 0.3, 
and frequency u) = 0.057 a.u. (800 nm). The gray curve rep- 
resents the HHG spectrum for a travehng wave. 
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FIG. 4: The same, as in Fig. 3, but for the frequency ui = 
0.043 a.u. (1054 nm) and i = 0.35. 
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FIG. 5: Harmonic emission rate for an ion with Za ~ 4, 
standing wave intensity 3.4 x 10^^ W/cm^ (^ = 0.4), and 
frequency u — 0.057 a.u. (800 nm). The gray curve represents 
the HHG spectrum for a traveling wave. 



Figures 3 and 4 display harmonic emission rate via 
modulus squared of the Fourier transform of the current 
density \Jn\^ as a function of the harmonic order N for 
an ion with Za = 3. The frequencies and relativistic pa- 
rameters of intensity are: uj = 0.057 a.u. = 0.3) and 
ui = 0.043 a.u. = 0.35), respectively. Figure 5 dis- 
plays harmonic emission rate for an ion with = 4 and 
field intensity 3.4 x 10^'' W/cm^ = 0.4) with the wave 
frequency: w = 0.057 a.u.. For the comparison we have 
also brought the spectra for a traveling wave which was 
obtained using the results of Ref. [18[. As we see from 



these figures, with the increase of the laser intensity the 
HHG rate with a standing wave field by many orders 
of magnitude is larger than HHG rate with a traveling 
wave. Besides, the relativistic quantum mechanical har- 



monic cutoff is shifted to the lower values compared with 
nonrelativistic one. 

Figures 6, 7, and 8 display harmonic emission rates 
when the relativistic parameter of the wave intensity is: 
^ = 1. In this case the relativistic effects are essential. 
Figure 6 corresponds to w = 0.184 a.u. and Za = 7. Fig- 
ures 7 and 8 correspond to Za = 5 with uj = 0.057 a.u. 
and u! — 0.043 a.u., respectively. For the last three figures 
the cutoff position is ~ 2.61Up, which is in good agree- 
ment with semiclassical analysis (see Fig. 2). For these 
parameters the harmonic emission rates with a traveling 
wave are negligible, meanwhile the relativistic rates with 
a standing wave field at photons energies approaching to 
MeV region, are considerable. 
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FIG. 6: Harmonic emission rate as a function of the harmonic 
order for an ion with Za ~ 7, = 1, and oj — 0.184 (248 nm). 
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FIG. 7: Harmonic emission rate as a function of the harmonic 
order for an ion with Za ~ 5, ^ = 1, and oj = 0.057 a.u.. 
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IV. CONCLUSION 



We have presented a theoretical treatment of the HHG 
by hydrogenhke ions with a large nuclear charge in the 
relativistic regime with a standing wave configuration, 
formed by the two linearly polarized counterpropagat- 
ing laser beams of relativistic intensities, at the distances 
much less than a laser wavelength. The investigation of 
the HHG with the such pump field configuration is con- 
ditioned by known fact that in the relativistic regime the 
significant HHG suppression takes place in the field of a 
strong traveling wave because of the magnetic drift of a 
photoelectron along the wave propagation direction. We 
have proposed a new scheme for relativistic HHG with 
an ultraintense traveling laser pulse and copropagating 
ions in a plasma. In this case the relativistic HHG oc- 
curs exactly without the magnetic drift of the photoelec- 
tron, because the traveling wave in the ions own frame 
of reference (moving with the laser pulse group veloc- 
ity) in underdense plasma is transformed into the purely 
uniform periodic electric field (in contrast to a standing 
wave - counterpropagating laser beams). Hence, we ex- 
pect that the efficient HHG with the current ultraintense 
laser pulses and forthcoming relativistic ion beams (at 
GSTDarmstadt accelerator) will be feasible in the plas- 
mas of conventional densities. 

On the base of the Dirac equation we have presented 
the analytic and fully relativistic quantum theory of 
HHG. With the help of the dynamic solution of the Dirac 
equation and Fast Fourier Transform algorithm we have 
calculated the harmonic spectrum. The obtained results 
have been applied to hydrogenlike ions with a moder- 
ate nuclear charge. The harmonic spectrum displays the 
significant difference compared both to non-relativistic 
and relativistic spectra with a traveling laser pulse. In 
particular, the shift of the cutoff position to the lower 
values of the harmonic order is considerable. In the rela- 
tivistic regime of interaction the harmonic emission rates 
for a traveling laser pulse are negligible, meanwhile the 
relativistic rates with a standing wave field, at photons 
energies approaching to MeV region, are considerable. 
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